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Abstract 

In the context of Type IIB LARGE volume orientifold setup equipped 
with poly-instanton corrections, the standard single-field poly-instanton 
inflation driven by a 'Wilson' divisor volume modulus is generalized by the 
inclusion of respective axion modulus. This two-field dynamics results in a 
"Roulette" type inflation with the presence of several inflationary trajecto- 
ries which could produce 50 (or more) e-foldings. The evolution of various 
trajectories along with physical observables are studied. The possibility of 
generating primordial non-Gaussianities in the slow-roll as well as in the 
beyond slow-roll region is investigated. We find that although the non- 
linearity parameters are quite small during the slow-roll regime, the same 
are significantly enhanced in the beyond slow-roll regime investigated up 
to the end of inflation. 
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1 Introduction 



On the way of understanding the early universe cosmology, the inflationary mech- 
anism has been proven to be quite fascinating as it successfully addresses several 
outstanding issues of the standard big bang scenario, e.g. the horizon problem, 
the flatness problem and the monopole problem. Initially, the idea of inflation was 
introduced to explain the homogeneous and isotropic nature of the universe at 
large scale structure [HE], however its best advantage is being utilized in studying 
the inhomogeneities and anisotropies of the universe, which is a consequence of 
the vacuum fluctuations of the inflaton (as well as the metric fluctuations). In this 
regard, the idea of inflation also provides a way to understand the physics which 
could be responsible for generating the correct amount of primordial density 
perturbations initiating the structure formation of the universe and the cosmic 
microwave background (CMB) anisotropies. The simplest (single-field) inflation- 
ary process can be understood via a (single) scalar field slowly rolling towards its 
in a nearly flat potential. The vacuum fluctuations of the inflaton result in an 
almost scale invariant spectrum with a small tilt reflecting unique predictions via 
the CMB radiation. Although the single-field inflationary models fit well with 
the current observation constraints [21 H], the present observational data is not 
sufficient to discriminate among the various models. In this regard, the detection 
of non-Gaussianity can be a crucial data to distinguish the various models in the 
ongoing/future experiments such as PLANCK [H El [TJ [8] . 

In the context of string cosmology, inflationary model building have been 
started quite early in [9]. Since string framework can provide several flat-directions 
(moduli), it is promising for the embedding of inflationary scenarios in string the- 
ory. The main requirement for this purpose is to identify a scalar which could 
play the role of the inflaton field, i.e. its effective scalar potential has to admit 
a slow-roll region. In this respect, those "moduli" that have a flat potential at 
leading order and only by a sub-leading effect receive their dominant contribu- 
tion are of interest. The perturbative effects [101 (TTJ as well as the instanton 
effects [121 [13] are proven to be extremely crucial,F especially for moduli stabi- 
lization purpose. As far as cosmological model building is concerned, the string 
inspired models could be given the hallmark of being '(semi)realistic' only after 
all the moduli could be stabilized and de-Sitter solutions be realized in [2]. Since 
then several dS realizing mechanisms have been developed in type IIB framework 
[T5l fT6l [TTf IT8l IT9l 120] . With the present understanding, it is fair to say that 
the moduli stabilization is quite well (and relatively much better) understood 
in type IIB orientifold models with the mechanisms like KKLT [TJ], Racetrack 
[2~T] [22] and the LARGE volume scenarios [23] • Significant amount of progress 
has been made in building up inflationary models in type IIB orientifold setups 
with the inflaton field identified as an open string modulus [2H [251 [261 [27] or a 
closed string modulus [281 ESI EH EU E2] . Along the lines of moduli getting lifted 
by sub-dominant contributions, recently so-called poly-instanton corrections be- 
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came of interest. These are sub- leading non-perturbative contributions which can 
be briefly described as instanton corrections to instanton actions. In the recent 
work [53], we have clarified the zero mode conditions for an Euclidean D3-brane 
instanton, wrapping a divisor of the threefold, to generate such a poly-instanton 
effect. Utilizing the poly-instanton effects, moduli stabilization and inflation have 
been studied in a series of papers [321 EH EH [36] . Meanwhile, the studies made in 
the context of axionic inflationary models in the type IIB orientifold framework 
have been quite promising too [211 133 ESI ESI SD1 El S21 S3] • In [S], inflation 
has been realized by a combination of the brane-motion and an axion while in 
[ISlHniHT],, it has been driven by a combination of a (divisor) volume mode and 
an axion. 

The signatures of non-Gaussianities are encoded in non-linearity parameters 
f NL (and r NL ) which parametrize the bispectrum (and trispectrum) and could be 
the direct evidence for detecting the non-Gaussianities. In fact, an observation 
of order one values for /jvl would rules out (most of) the single field inflationary 
models H. and thus can be a potential discriminator towards picking up a more 
promising one. This boosts up the motivation for the study of non-Gaussianties 
in recent years, although it has been initiated very early |52j . A review on initial 
attempts regarding /jvl computations along with observational constraints can 
be found in [5]. In the standard single-field models, the /nl parameter is usually 
suppressed by slow- roll parameters (e, 77), e.g. Jnl = n s — l has been reported in 
|52j . Subsequently, the multi-field inflationary models started getting significant 
attention but large /tvl could not be realize in the initial attempts [531 [Ml ESI EE] • 
For generating observable non-Gaussianities, enormous amount of work has been 
done in recent couple of years and several explicit multi-field models with large 
/jvx values have been found. These multi-field models can be categorized to be 
of separable and non-separable types. 

For the separable type models, either the inflationary potential [531 Ell EHl 
I5T1 EB] or the Hubble rate [521 ED] is of sum or product separable form, and 
several models are available now which can produce large non-Gaussianities @. 
The possibilities of generating detectable large values of Jnl during the evolution 
of two- fields within the slow- roll region has been proposed in [631 EH E5]. In 
these models, non-Gaussianities are generated after horizon exit and consist of 
non-adiabatic perturbation modes. It implies such non-Gaussianities to be of 
local-type and hence distinguishable from the other shapes of non-Gaussianity 
which are realized during horizon exit [661 EZ] ■ 

Unlike the separable type models, the non-Gaussianities issues in the non- 
separable type of potential are relatively less studied. In [561 EH ESj, a con- 
cise analytic formula for computing the non-linear parameter for a given generic 

3 There has been some proposals of single-field models with large Jnl values [48l l49l l50l loT) . 
A nice recent review about generating large non-Gaussianities for separable type setups 
can be found in [STJ [52] . 
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multi-field potential has been developed for both the slow-roll region and beyond 
slow-roll region. The possibility of generating large non-Gaussianity parameter 
has been shown by 'momentarily' violating the slow-roll conditions. Following 
the ideas developed in [68J , order one values of the f^i parameter has been re- 
ported in the beyond slow-roll region for a two-field large volume axionic setup 
[70] . Recently, some examples with (non-) separable multifield potentials have 
been studied in [71] which can produce large detectable values for the non-linear 
parameter f^L and t^l- Higher order correction to the non-linearity parameters 
have been studied in |7J1 Efl EH E5] . 

In the framework of type IIB orientifolds, several single/multi-field models 
have been studied for aspects of non-Gaussianities. The curvaton scenarios have 
been proposed in [76] for the Kahler moduli inflation setup in the context of 
LARGE volume scenarios. A case study have been performed in [77] for a blow- 
up multi-field inflationary setup (resulting in f NL ~ 0.01) and a curvaton scenario 
(resulting in large f^L values). Later on, detectable large values of f^L have also 
been proposed via the modulated reheating mechanism in [36]. The computation 
of non-Gaussianties in racetrack models has been made in [78] and within the 
framework of KKLT-like setup, a two-field inflationary model has been proposed 
with inflaton dynamics governed by the Calabi Yau volume mode and the respec- 
tive C4 axion which complexifies the divisor volume mode [15]. This idea has 
been extended in the context of LARGE volume scenarios by generalizing the 
Kahler moduli (blow-up) inflation with inclusion of the corresponding C4 axion 
in [161 H7J. The presence of various types of trajectories resulting in 60 (or more) 
number of e-foldings have been reported in this so-called "Roulette inflation" 
model and a subsequent investigation of non-Gaussianities in such a setup has 
been made in [79] resulting in small values of the fNL parameter. 

In this article, we present a systematic analysis for a two-field inflationary 
model driven by a combination of a Wilson divisor volume modulus and an axion 
modulus in the Poly-instanton setup. The so-called Wilson divisor has a single 
non-trivial Wilson line modulino and is relevant for generating poly-instanton 
contributions which will be summarized in Section |2j For analyzing the non- 
Gaussianities, we proceed with the formalism developed in [681 EH] which is valid 
for a given generic multi-field potential in beyond slow-roll region also. We find 
that although the non-linearity parameters are quite small during the slow-roll 
regime, the same are significantly enhanced in the beyond slow-roll regime. The 
two-field dynamics is such that the non-linearity parameters (/nl^nl and gNh) 
do not get frozen at the horizon exit and keep evolving up to the end of inflation 
where it acquires a large value. Here, we stress that following the motion/decay 
of inflaton fields after the end of inflation and addressing reheating issues etc. can 
be extremely crucial in such models. However, the same is beyond the scope of 
the present analysis, and for the time being we assume that the respective values 
for these non-linearity parameters realized at the end of inflation do not change 
significantly before and by the reheating process. 
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The article is organized as follows. In Section [21 we start with a brief review 
of the relevant background about the poly-instanton setup [33] and the related 
moduli stabilization part [32]- In Section [31 we generalize the single field poly- 
instanton inflationary process into a two-field inflationary process with the inclu- 
sion of the corresponding C4 axion and discuss the subsequent relevant changes in 
the slow-roll parameters. In Section HI we present a detailed and systematic (nu- 
merical) analysis of evolutions of various inflationary trajectories (as well as the 
other physical observables) in terms of number of e-foldings. This analysis shows 
that there is a "Roulette" type inflation. In section \5\ we investigate the pos- 
sibility of generating finite/detectable values for the primordial non-Gaussianity 
parameters /tvl? Ttvx; Qnl and observe that although in the slow-roll region these 
are negligibly small, beyond it there exists the a possibility of enhancement to 
large values. Finally, in section E] we give our conclusions followed by an appendix 
providing some of the lengthy intermediate expressions. 

2 Poly-instanton Setup 

In this section, we collect the relevant ingredients for generating the poly-instanton 
corrections in a type IIB orientifold setup developed in [33] and briefly summarize 
the moduli stabilization mechanism discussed in [32]. Building on the same, we 
continue with the study of a generalized two-field inflationary process as well as 
with the computation of the non-linearity parameters (e.g. fjsiLi t nl and gNLj- 

2.1 Poly-instanton corrections 

Now, we recall some results from [331 180] (see also [81] for related studies) on 
the contribution of poly-instantons to the superpotential in the framework of 
Type IIB orientifold compactifications on Calabi-Yau threefolds with 07- and 
03-planes. In this case the orientifold action is given by Qa(— 1) Fl , where a is a 
holomorphic, isometric involution acting on the Calabi-Yau threefold Ai. 

The notion of poly-instantons [3^1 180] means the correction of an Euclidean 
D-brane instanton action by other D-brane instantons. The configuration we 
considered has two instantons a and b with proper zero modes to generate a 
non-perturbative contribution to the superpotential of the form 

W = A a exp~ Sa + A a A b exp~ Sa ~ Sb + ... , (1) 

where A a ^ are moduli dependent one-loop determinants and S a} b denote the clas- 
sical D-brane instanton actions. 

In Type IIB orientifolds models, the sufficient conditions for the zero mode 
structures for poly-instanton corrections to the superpotential have been worked 
out in [33] (for orientifold with 05- and 09- plane, see [80J). Both instantons, 
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a and b, should be 0(1) instantons, i.e. a single instanton placed in an ori- 
entifold invariant position with an 0(1) projection which corresponds to an 
SP— type projection for a corresponding space-time filling D7-brane. Instan- 
ton a is an Euclidean £"3 instanton wrapping a rigid divisor E in the Calabi-Yau 
threefold with H lfi (E,0) = H 2 ' (E,O) = while instanton b is an Euclidean 
i?3-brane instanton wrapping a divisor which admits a single complex Wilson 
line Goldstino, i.e. a so-called Wilson line divisor with equivariant cohomology 
H*'°(W, O) = (1 + , 1 + , 0) under the involution a. In fact, the sufficient condition 
for a geometric configuration to support the poly-instanton correction is that it 
precisely contains one Wilson line modulino in H+(E, O). 

Some concrete Calabi-Yau threefolds both with and without K3 fibration 
structure were presented in [33] which featured all the requirements mentioned 
above. For the present purpose, we focus on one particular model which not only 
admits a Wilson line divisor W but also rigid and shrinkable del Pezzo divisors. 
Such geometries are particularly interesting for studying moduli stabilization as 
they give rise to a swiss-cheese type Kahler potential. The Calabi-Yau threefold 
Ai is given by a hypersurface in a toric variety with defining data, 
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with Hodge numbers (h 21 ,h u ) = (72,4) and the corresponding Stanley-Reisner 
ideal 

SR = {x 1 x 2 , x 4 x 7 , x 5 x 7 , xt x 4 x 8 , x 2 x 5 x 6 , x 3 x 4 x 8 , x 3 x 5 x 6 , x 3 x 6 x$} . (2) 

This geometry admits two inequivalent orientifold projections a : {x 4 -B- — £4, x 7 -R- 
— £7} with h}}(A4) = so that for the Wilson line divisor W = Dg = {xs = 0} 
the Wilson line Goldstino is in H+(W, O). It was checked that the D3- and D7- 
brane tadpoles can be canceled. For the analysis in the following sections, we 
focus on the involution x 7 -H- — x 7 . The corresponding topological data of the 
relevant divisors are shown in table [TJ 

We also showed that there are no extra vector-like zero modes on the intersection 
of E3(lD7, i.e. all sufficient conditions were satisfied for the divisor W to generate 
a poly-instanton correction to the non-perturbative superpotential 

W = A x exp (-2ttTi) + A x A 8 exp (-2tiTi - 2ttT 8 ) + 

A 7 exp (-a 7 T 7 ) + A 7 A s exp {-a 7 T 7 - 2nT s ) + . . . . 

Taking into account the Kahler cone constraints, the volume form for this model 
could be written in the strong swiss-cheese like form 

V = § (j-( n + 3r 6 + 6r 7 + 3r 8 ) 3 / 2 - V2rf /2 - 3r? /2 - 3(r 7 + r 8 ) 3 / 2 ) . (4) 
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divisor 


{h™,h w ,h w ,h 11 ) 


intersection curve 


D 7 = dP 7 


(l+,0,0,8+) 


W : C g=1 




(1 + ,0,1+,21+) 


W : C g=l 


D 8 = W 


(l+,l+,0,2 + ) 


D 5 : C g= i, D 7 : C g= \ 


Dj = P 2 


(1+, 0,0,1+) 


D$ : C g= o 



Table 1: Divisors and their equivariant cohomology under x 7 <-> — x 7 . The first two 
lines are 07-plane components and the remaining two divisors can support E3 instan- 
tons. The D 7 divisor also supports gaugino condensation. 

This volume form indicates that the large volume limit is given by tq — > oo while 
keeping the other shrinkable del Pezzo four-cycles volumes r 17 and the Wilson 
line four-cycle volume r 8 small. In fact, all the models studied in [53] shared a 
similar strong swiss-cheese like volume form with the same intriguing appearance 
of the Wilson line Kahler modulus. 



2.2 Moduli stabilization 

A generic orientifold compactification of Type IIB string theory leads to an effec- 
tive four- dimensional J\f = 1 supergravity theory^. In the closed string sector, the 
bosonic part of the massless chiral superfields arises from the dilaton, the com- 
plex structure and Kahler moduli and the dimensional reduction of the NS-NS 
and R-R p-form fields. The bosonic field content is given by 



T = C {0) + ie~^ , IP = v? + iv\ i = l...h 



21 

+ ' 



G a = c a -rb\ a = l,...,h 1 }, 

1 i 
T a = -Ka^f + i(p a - K-aab c"& 6 ) + - TK aab b a b b and a = 1, . . . , h+ , 



(5) 



where c a and b a are defined as integrals of the axionic and forms and 
p a as integrals of over a basis of four-cycles D a . From now on, as in our 
aforementioned concrete example, we assume h]} = 0. 

The supergravity action is specified by the Kahler potential, the holomor- 
phic superpotential W and the holomorphic gauge kinetic function. The Kahler 
potential for the supergravity action is given as, 



K = -rn( -i(r-f) 1 -In I -i J ft A fij - 2 \n[V(T a )j , (6) 

where V = ^K a ^r f t a t^t y is the volume of the internal Calabi-Yau threefold. The 
general form of the superpotential W is given as 

W= I G 3 An + y2A E (r,U i )e- aE ^ Ta (7) 

JM TP 



3 For a review on moduli stabilization and relevant compactification geometries, see [821l83j . 
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with the instantonic divisor given by E = Yll a D a - The first term is the 
Gukov-Vafa-Witten (GVW) flux induced superpotential [TTJ and the second one 
denotes the non-perturbative correction coming from Euclidean D3-brane in- 
stantons (ag = 2tt) and gaugino condensation on U(N) stacks of Z)7-branes 
(oe = 2tt/N) [12]. In terms of the Kahler potential and the superpotential the 
scalar potential is given by 



V = e K \^K iJ V I WV 3 W -?>\W\ A j , (8) 

where the sum runs over all moduli. As usual in the LARGE volume scenario 
|23j . the complex structure moduli and the axio-dilaton are stabilized at order 
1/V 2 by K and the GVW-superpotential respectively. Since the stabilization of 
the Kahler moduli is by sub- leading terms in the V -1 expansion, for this purpose 
the complex structure moduli and the dilaton can be treated as constants. 

Let us proceed with the following ansatz for Kahler and superpotential) which 
is well motivated by our mathematical discussion in the previous subsection 12. 1^ 

K = -2 my, 

W = W + A s e~ asTs + A S A W e ~ a ^- a ^ (9) 

JD „—bsT a D R „—bsT a —b w T w 

u s e u s u w e , 

where y = V{T a ) + C a > such that 

3 

y = i b {T b + f b )l- UTs + T s )l - U((T S + T.) + (T w + T w ^j 5 + C a , . (10) 

Note that, for h 1 } = 0, the M = 1 Kahler coordinates are simply given as 
T a = T a + ipa- Here, C a > denotes the perturbative a' 3 -correction given as [10] 



K iJ VjWVjW -3\W\ 2 ) , 



C a , = _x(M) { r-f)ip) (n] 



4(2vr) 3 (2z)5 



with x{M) being the Euler characteristic of the Calabi-Yau. The large volume 
limit is defined by taking r& — > oo while keeping the other divisor volumes small. 

In the large volume limit, the most dominant contributions to the generic 
scalar potential V(r&, r s , r w , p s , p w ) are collected by three types of terms. In the 
absence of poly-instanton effects, they simplify to 

V(V, r s , p 8 ) ~ V Q ,(V) + V npl (V, r s , p.) + V np2 (V, r a , p s ) , (12) 



6 Wc consider the racetrack form of superpotential as it has been realized in [32] that the same 
is needed for having a (non-susy) minimum which could be trusted in the effective supcrgravity. 
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where 
V* = 

v npl = 

V np 2 - 



3C a ,\W \ 2 
2 V 3 ' 

4 a s A s W e~ a ° Ts r s cos( a s p s ) 4 b s B s W e~ bsTs r s cos( b s p 



(13) 



V 2 V 2 

2V2a 2 s A 2 s e- 2a ^^/T s 4^ a s b s A s B s e ~K+M^ v ^cos(( a s - b s )p s ) 



+ 



2V2b 2 B 2 e- 2b ° T °JT s 



36V 



As expected, the potential ( lT3"j) does not depend on the Wilson line divisor volume 
modulus r w . Therefore, at this stage it remains a flat-direction to be lifted via 
sub-dominant poly-instanton effects. 

The extremality conditions dy\ = d Ts V = d Ps \ = are collectively given as 



Wo 



a s Ps 



V {b s X 2 - a s Xi) 


b s X 2 {-l + 4:b s r s ) - a s Ai(-l + 4a s r s ) 




M 2 (-l + b s T s ) - a s Ai(-l + a s r s ) 





32V2£ s T${b 2 s \ 2 -a 2 s \ 1 ) b s \ 2 (-l + b s r s ) - asXii-l + a s r s 



(14) 



a s Ai(-l + Aa s r s ) - b s X 2 (-l + Ab s r : 
Ntc with Ai = A s e" asTs and A 2 = B s e~ b ^ ! 



I 2 



One finds that t s gets stabilized in terms of C a i as r s ~ (C a i)s and then V gets 
stabilized via an exponential term exp(a s r s ) (encoded in Aj's) so that the overall 
volume of the Calabi-Yau threefold is exponentially large. 

Now, in addition to the leading order standard racetrack terms (TIB"]) , the 
generic scalar potential also involves sub-dominant contributions, which are fur- 
ther suppressed by powers of exp(— a w r w ). Collecting these sub-leading terms, 
the effective potential for r w ,p w becomes 

V(V, T s , T W , Ps , Pw ) = V LVS (V, T t , p s ) + V(T W , p w ) , (15) 

where, V LV (V,t s , p s ) is the racetrack version of the standard large volume po- 
tential which stabilizes the Kahler moduli t\> (or V) and t s at order V -3 and 

4W 



'V\Tyj) Py 



V 2 



+ 



Ai A w e amTw (a s Ts + a w T w ) cos(a w p w ) 

- X 2 B w e ~ bwTw (b s T s + b w r w ) cos(b w p 
4V2(6 a A 2 -a a Ai) A /7i 



36V 



(16) 



X 2 (b s - b w )B w e bwTw cos(b w p w ) 
- Ai(a a - a w )A w e~ awTw cos(a w p w ) 
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. After stabilizing the heavier moduli V, r s , and p s -axion at their respective 
minimum and using a w = b w along with eliminating Wq via the first relation in 
eq. (j!4p . the above effective scalar potential can be written as 



V(t w , pw) = V + e awTw (ji! + fa t w ) cos(a w p % 



(17) 



Here V , pi, p 2 are constants depending on the stabilized values of the heavier 
moduli as 



A 4 ! — Po 



4r.(( 



a, — a, 



; )A w \i — (b s — a w )B w X 2 ) 



+ 



p 2 — 



r s (b s B w X 2 - a s A w \i) (a s Ai(-l + Aa s r s ) - b s X 2 (-l + 46 s r s )) ' 
a s Xi (-1 + a s T s ) - b s X 2 (-l + b s T s ) 

(B W J 2 - A^Ai) (a s Ai(-l + 4a s r s ) - 6 S A 2 (-1 + 46 s r s )) 



a s Ai(-l + a s r s ) - b s \ 2 (-l + b s r s ) 



with 



Ik) 



\/2( a s Xi - b s X 2l 



(19) 



After stabilizing the p^-axion at its minimum ~p w = 0, it has been shown (in 
that the divisor volume modulus r w corresponding to the Wilson divisor D w gets 
stabilized to 



1 ^ 

1 - CLw 

M2 



(20) 



where ^ < is required for stabilizing r w inside the Kahler cone. A couple 
of benchmark models have been constructed for a set of sampling parameters 
displacing the Kahler modulus t w away from its minimum, investigations of in- 
flationary behavior have been made in a single field approximation [32]. In this 
article, we generalize the inflationary model via including the dynamics of corre- 
sponding p w axion and discuss various subsequent cosmological implications. 



3 Two-Field Poly-Instanton Inflation 

In all the remaining sections, we will be assuming that the heavier moduli are 
stabilized at their respective minimum position and we consider the effective po- 
tential for the two lighter fields (t w and p w ) given by (fl7|) . Furthermore, we 
assume that a suitable uplifting of the AdS minimum to a dS minimum can 
be processed via an appropriate mechanism (e.g. by the introduction of anti-D3 
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brane |14) 7 I or other mechanism like using dilaton-dependent term in superpoten- 
tial [19|). Thus, the resulting inflationary potential looks like 

V in{ {r w , p w ) = V up + V + e~ awTw (//i + /i 2 r w ) cos(a wPw ) (21) 

Now, we state the effective potential studied before in roulette inflation [IB] which 
after stabilizing all but one Kahler moduli effectively looks as, 

if , , v , v , V^e~ 2anTn , r n e- a - T - cos(a ra p n ) 

VrouletteXJn, p n ) — Vup + Vq H = 1 =5 

where r n is stabilized at r n ~ InV. Let us mention some important differences 
between the aforementioned potentials. First, the geometric origin of these two 
potential are different as we have explained in the previous sections. The former 
one comes from a Wilson divisor volume modulus while the later one comes from 
a blow-up volume modulus. Second, these moduli are stabilized by corrections 
which scales differently in terms of CY volume. Moreover, the stabilized values 
at the respective minimum have different volume scaling; r n is stabilized at order 
InV while r w is stabilized by eq. (l20|) . Here, the V dependence, which could have 
appeared via r s , gets effectively canceled. 

The uplifting term V up in eq. ([2~T]) needs to be such that the uplifted scalar 
potential acquires a small positive value (to be matched with the cosmological 
constant) when all the moduli sit at their respective minimum. This potential 
has the following set of critical points 

(I) . t w = , a w p w = 2pn (23) 

(II) . r w = , a w p w = {2p+ 1)tt 

a w P2 

where peZ. Moreover, in order to trust the effective field theory we need ^ < 0. 
In order to ensure the minimum, one has to consider the Hessian V a b (evaluated 
at these two sets of critical points) which is given as 

(l I "TO Ml \ 
Ta w p 2 e "2 \ , s 

1 | a w Ml , 
Ta w fi 2 e "2 ) 

where =F sign corresponds to critical points {%) and (ii) respectively. Hence de- 
pending on whether {jii < 0,p 2 > 0} or {jix > 0,p 2 < 0}, we find that one 
set of critical point corresponds to minima while the other corresponds to max- 
ima. From now on, we fix our notation with a sampling of parameters such that 

7 For uplifting process using anti-D3 brane, there has been some sensitive issues as mentioned 
in [H|85l[86]. 
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> 0,yU 2 < 0} (as in [32]) and performing the redefinitions r w = (pi, p h 
the uplifted scalar potential takes the form 



infill; V2J 



9s 



-i+- 



+ e 



(A*i + A*2 0i) cos(a w 2 



.(25) 



Here, a proper normalization factor (|^) e Xcs has been included [28}, where i^cs 
denotes the Kahler potential for the complex structure moduli. We assume that 
e K cs ^ (9(1). Furthermore, we set the numerical parameters for moduli sta- 
bilization similar to the ones chosen in one of the benchmark models (in [52]). 
The parameters which would be directly relevant for further computations in this 
article are, 



= 2.9 x 10" 8 , Hz = -1.9 x 10~ 8 , a w = 2n,g s = 0.12 (26) 
V = 905,r a = 5.7,e« o = l/(6V2) 

The 'effective' non-flat moduli space metric relevant for inflaton dynamics 
can be computed as 



36 



lab 



2V2V VTs+<Pl 







3g a 



(27) 



2 V2VVt s +0i 



Utilizing the large volume limit, the same is written from the Kahler metric 
component K T 5? expressed in the real moduli basis {0i,</>2}. The Christoffel 
connections are defined as T a bc = Q ad T bc d = \Q ad ipbQcd + dcQbd — ddGbc) and the 
only non-zero connection components are 

1 



r 1 
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4(r 3 + 0i) 
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r 2 r 1 

1 21) 1 22 



4(r 3 + 



Furthermore, the non-zero components of the Riemann tensor, defined as R 

Uc*- db u d L cb* 1 ck L db 1 dk L cbi dIt; 



bed 



R\ 
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2(r 3 



p2 pi 
- rL 121' ' 221 



2(r 3 
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Under the sampling (T2"6"l) . the effective two-field inflationary potential from eq. ( 12 1 p 
is shown in FigHJ 



4 Evolution of Trajectories 

In this section, we investigate the field evolutions and look for the possible infla- 
tionary trajectories which could reproduce 50 (or more) e-foldings. For a given 
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Figure 1: The effective potential as a function of the moduli t w and p w . 

initial condition, we numerically trace the entire inflationary trajectories includ- 
ing beyond slow-roll region. Using the background N e-folding number as the 
time coordinate, i.e. dN = Hdt, the Einstein-Friedmann equations are obtained 

as 

d 2 „ n d6 b d(b c ( 1 dH\ d6 a G ab d b V , , 



dN 2 ^ dN dN V H dN J dN H 2 



Using expressions f)28al) and (I28bl) . one can derive another useful expression for 
variation of Hubble rate in terms of e-folding, 

ldH V -3. (29) 



HdN H 2 

For numerical convenience, we solve these equations in the time basis t and then 
change the result back to the basis N e-folding. As introduced in we will 
follow the field redefinitions given aa^| 



^ = a , ^ = 0- =(_?_), where a = 1,2. (30) 



8 The use of this notation would be more clear in the next section regarding computations 
of non-linearity parameters. Further, we will be using a combined indexing A such that any 
object O a has two components given as C" 4 = {Of, 0%}. 
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which translates the second-order background equations of motions eq. f)28ap into 
two first-order ODEs as follows 

Ft - ** 





1 


\dN ) 




-3^ 


r ab V b 

~ Q H 



dN 

r?a _ Dip 2 

= ~dN~ = ^ ~ Q H- (31) 

where D is the covariant derivative defined as Dip?, = dip?, + T a i lc <p\d(p\ , subject 
to the constraints H 2 = | [V + \ G a b^2 y^) • Then eq. (l29|) will be simplified as 
^ = — \ Q a \ 1 ip? 1 ip\. Now, in the context of studying inflationary aspects, one 
has to look at the sufficient conditions for realizing slow-roll inflation which are 
encoded in the so-called slow-roll parameters. For multi-field inflationary process 
with inflatons moving in a non-flat background, these slow-roll parameters are 
given as, 

1 dH 1 de 

e = -Jp^ ^THdf (32) 

Now, we can solve the background field equations (13TT) to get the full trajec- 
tories under different initial conditions. We choose a (O) = 0q and ^p^p\t=o — 
0; for a G {1, 2} as a set of initial conditions and trace the corresponding trajecto- 
ries up to the end of inflation (see Table [2]). Figure [2] shows the complex evolution 
of trajectories for some samples of general starting conditions where the value of 
iV e-folding at the end of inflation is labeled on each of these trajectories. Note, 
although it looks like a single field inflation in the late period of trajectory, the 
axion field in fact oscillates (Figure [3]) during that time before settling into its 
minimum. The fact that both of the fields in each trajectories keep evolving up 
to the end of inflation might be relevant for generating large Non-Gaussianity. 
This indicates that isocurvature modes are not completely finished even though 
it looks as a single field motion in the final stage of the trajectories. The iV 
e-folding when the axion p w first crosses its minimum is denoted as N' . One has 
to note that at this time TV', the Wilson divisor volume mode is still quite far 
from its respective minimum. However, it rolls quickly enough to reach there. 
All of these points are already in the beyond slow-roll region as the slow-roll 
parameter 77 > 1. 

The various inflationary trajectories shown in Figurej2]can be classified in the 
following categories 

(a) If the axion initial condition is such that the axion is minimized at its respec- 
tive minimum, then two-field inflationary process reduces to its single field 
analogue which has been studied in [52] • These are stable trajectories and are 
attracted towards the respective valley in a straight line like the trajectory in 
Figuref2]with Np = 62. These can produce the required number of e-foldings 
if the Wilson divisor volume mode is displaced significantly away from the 
minimum. 
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3.5 4.0 4.5 5.0 5.5 



Figure 2: The full inflationary trajectories for various initial conditions where the 
value of N e-folding at the end of inflation is labeled on each of these trajectories. 
Various minima in dark blue are separated by maxima in light blue shade. 

(b) If the axion initial condition is a little bit away from the minimum, the tra- 
jectories rolls to the nearest valley and trace towards the respective minimum 
like those trajectories in Figuref2] with Np = 1, 53.4, 65, 98.8. 

(c) If the axion initial condition starts with its value at the maximum, this results 
in an unstable trajectory directed straightly outwards from the respective 
attractor point showing a run- away behavior like the yellow trajectory in 
Figure [21 

(d) The trajectories starts from axion initial conditions being closer (but not 
exactly equal) to some maximum value as well as the initial values for the 
divisor volume mode being not very far from its respective minimum, one 
observes that inflationary trajectories cross several axion-ridges before getting 
attracted into a valley. This can be understood from the fact that this class of 
initial conditions is such that the initial potential energy is just a little higher 
to begin with (as shown in FigureJH) and the N e-folding increase very slow 
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Figure 3: Oscillation of axion-field p w during late period in the beyond slow- 
roll regime. The N e-folding when the p w crosses its minimum is denoted as 
A^t = {64.67, NULL, 53.08, 98.42}. At the time the divisor volume mode 
is still quite far from its minimum and thus both fields are dynamical and it is 
already in beyond slow-roll regime. Notice that the second trajectory corresponds 
to a single-field inflation. 

at the beginning of these trajectories, see FigurefS] with N F = 7.4, 44.2, 270. 




Figure 4: The potential V as a function of e-folding N, where the value of AT 
e-folding at the end of inflation is labeled on each of these trajectories. 

For most of these trajectories except the single-field one, they momentarily 
undergo some sort of quick-roll region before staring the slow-roll (see Figure 
HI [5] and [6]). The slow- roll region ends when any of the slow roll conditions 
e -C 1, rj <§C 1 is violated. The evolution of the slow- roll parameter e and r\ 
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Pw 


AT 
l\p 


a 


5 


1 


62 




4 


0.3 


1 


b 


4.5 


1.475 


53.4 




4 


0.496 


65 




3.9 


1.495 


98.8 


c 


3.5 


-0.5 






3.65 


0.2 


7.4 


d 


3.4 


0.3 


44.2 




3.7 


0.4 


270 



Table 2: Initial conditions for these trajectories shown in Figurej2j 

depending on the variation of e-folding for four of these trajectories are shown 
in Figure [5] from which one can see that these parameter changes dramatically 
at the end of inflation. Also, there is a region in field space where there is a 
strong violation of slow-roll condition via rj ^> 1 before the end of inflation. This 
beyond slow-roll region can be interesting from the non-Gaussianities point of 
view as explained later. Furthermore, we can define the scalar perturbation power 
spectrum P s (in the slow-roll region) as P s = / 2 ^2 e . Figure |6] shows the evolution 
of the Hubble rate H and power spectrum P s with e-folding variations for the 
corresponding trajectories. It shows that the Hubble rate H is almost constant 
(at 10 _8 M p ) during entire inflationary process including the beyond slow-roll 

region also. This indicates a high scale inflation as M in f ~ ~ H 2 ~ 10 14 GeV. 
The power spectrum for scalar perturbations P s reaches a value of order 10 -9 
at the horizon exit. Also, within the slow-roll limit, the spectral index is found 
to be negligibly small. All of these results are consistent with the observational 
constraints today. 

5 Primordial Non-Gaussianities 

The signatures of non-Gaussianities are encoded in a set of non-linearity param- 
eters which are commonly denoted as fNL, T NL and Qnl- These are generically 
related to the n-point correlators of curvature perturbations; the 2-point corre- 
lators simply give rise to a Gaussian shaped power spectrum while the 3-point 
correlators are related to the bi-spectrum which encodes the non-Gaussianities 
via the non-linearity parameter /jvl- Similarly, the 4-point correlators give rise 
to a tri-spectrum via tjvl and g^i parameters. These non-linearity parameters 
can be computed in the 5A-formalism which relates the curvature perturbations 
to the difference between the number of e-foldings 5N of two constant time- 
hypersurfaces [52] . 

C(t,x) ~5N = H5t 
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N 

Figure 5: Slow-roll parameter e and n as a function of e-folding N, where the 
value of N e-folding at the end of inflation is labeled on each of these trajectories. 
The enhancement of the rj parameter towards the end of inflation is illustrated 
for one trajectory corresponding to Np = 65 in the last plot. 

Following the redefinitions of the background field evolutions as made in the 
previous section, the perturbations of the scalar field on N = constant gauge can 
be expressed as, 

6(p A (X, N) = tp A (X + 5X, N) - <p A (X, N) 

where A's are In— 1 integration constants (for an n- component scalar field) which, 
along with N, parametrizes the initial values of the fields [68| [69] . The curvature 
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Figure 6: The Hubble rate H and the power spectrum P s as a function of e- 
folding N, where the value of N e-folding at the end of inflation is labeled on 
each of these trajectories. 



perturbations at each spatial point of the field space are subsequently expressed 
in terms of variations of the number of e-foldings in various field directions, 



N 



C(iV F ,x) = J]liV^ 2 ....^^(x)^(x)... 

d n N(N F , V A ) \ 



.Sip 



An 



(33) 



A\Ai....An 



where ip A ^ corresponds to an unperturbed trajectory and Np corresponds to a 
final time-hypersurface of uniform energy density. The non-linearity parameters 
In Li t~nl and g^L are generically defined as, 

5 N A N B N AB 



1 



NL 



T~NL 



9NL ~54 



6 (N V N V ) 2 ' 
N A N AB N BC N c 

25 N A N B N c N ABC 



(34) 



(N"N V ) 



where the field variations of N are defined as N A = d A N, N AB = d AB N and 
N A = Q Nb- Now, the main task in computing the non-linearity parameters 
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is to find out the fields variations of the number of e-foldings and for the same 
we would follow the strategy developed in EH] using 5N formalism. It is 
important to mention that this approach is valid in the beyond slow-roll regime 
as well. So one can explore the evolutions of the non-linearity parameter in the 
non slow-roll regime up to the end of inflation. 



5.1 Strategy for computing /nl^nl an d gNL 



Let us briefly discuss the necessary ingredients from [HSl EH] which are relevant for 
computing the non-linearity parameters. For a given generic form of the scalar 
potentiaL the concise expressions for the non-linearity parameters Jnl,tnl and 
9nl aref 



NL 



6 (iV* Q?y 



1V AB U F U F 



N F 



N A Qbc @ S ® C dN 



(35a) 



TNL 



A AB n A (N*)n B (N*) 



(35b) 



25 

9NL — TT 



54 {Nie?y 



N% EC eie B F e c F + 3 / n A Q$ c e B e c dN 

Np 



N A Qbcv © B © C ® V dN 



N, 



(35c) 



Let us explain the meaning of the various symbols, 



It is important to recall that in f l35a|35b|35cp . any object O a has two 
components O a = {Of, 0%} and the indices are appearing due to field 
redefinitions fl30l). i.e. 



dt 



where a G {1, 2}. 



• The expressions of various derivatives {N A (N), Q A (N), Q A (N)} at a generic 
time during the inflationary dynamics is required for computing the inte- 
grals. These vector quantities have to be computed by solving the respec- 
tive ODEs. This is among the main task of the computation and we will 
elaborate more on it later. 

9 To avoid any possible confusion with the notations, we mention that all capital letters 
appearing as the sub/super scripts are denoted in calligraphic font. 
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The symbols Q A (N) are defined as 

e- 4 (TV) = (N,N*) Nf ■ where (36) 



Texp^ P A B {N)dN 



with A A B (N*,N*) = 5£. 



In the above, T denotes time-ordering. The expressions for P-^(A r ) are 
given as, 



>ai 1 



n = + ^3 - ^ ^ a (37) 

The symbol A AB is defined as, 

(6<ft6tf)=A M (§) 2 . (38) 

In general, A AB depends on the non-flat background metric. The respective 
expressions are given in appendix (jAj) taking the slow-roll corrections [87] 
into account. 

The expressions of various derivatives of e- folding N evaluated at some final 
constant time-hypersurface ip (e.g. iVj[, N^ B , N^ BC ) are given as, 



N 



H 



A 



A ~ \H V F V 



W i = "(^l (39) 



at v =<p(°)(N F ) 



at ip= V (°)(N F ) 

F I ZABC 



N A bc 



The expressions for quantities H^(N), Hab(N), H^sc(N) ,U^b(N) ,Z^bc{^) 
as well as Q A isc(^) an d Q^bcvi-^) involve various derivatives of the scalar 
potential and the Hubble rate. Being quite lengthy, their explicit expres- 
sions can be found in (TA~1) . 



The main advantage of the formulation involving the redefinitions fl30|) is that 
this simplifies the computation of the non-linearity parameters /jvx> tnl and g^L- 



21 



It remains to solve first order ODEs for vector quantities (like Q A (N), N A (N), 
Q A (N)). Furthermore, this formulation reduces the number of 0(n 2 ) calculations 
to 0(n) where n is the number of scalar fields involved in the dynamics. This 
makes numerical calculations much more efficient in a multi-field scenario which 
has a large number of scalar fields. 

The expressions for 0- 4 (iV), N A (N), Q A (N) required for solving the integrals 
can be obtained by solving the following set of first order ODEs, 

^N A (N) = -P A B (N)N B (N), 

-^Q A (N)=Pi(N)Q B (N), (40) 

§jtl A (N) = -n B (N) P B A (N) - N B {N) Q B AC (N) Q C (N). 

Each set of these involves four ODEs for each quantity Q A (N), N A (N) and 
Q A (N). Of course, even this simplification of the problem into solving first or- 
der ODEs does not allow to proceed analytically for a given generic multi-field 
potential. However for concrete models, it is much easier to solve the afore- 
mentioned twelve ODEs numerically. As an important observation, the first two 
expressions of fj4*0|) are mutually dual to each other. This implies that the quan- 
tity X(N*) = N A (N*) ©^(iV*) is constant irrespective of N*. The algorithmic 
approach of solving the ODEs in (l4"0"j) can be summarized as, 

• First, one numerically solves the set of ODEs for N A (N) by using the 
boundary conditions corresponding to the final constant time-hypersurface 

Then, utilizing the numerical solutions 



N A (N F ) = M = 
obtained, one traces back to N A (N*). 



H.A 



Using the backward traced values for N A (N*), one gets the initial conditions 
O^-ZV*) = A AB (N i „) 7Vg(7V*) and subsequently one numerically solves the set 
of ODEs for Q A (N). 

Utilizing the solutions for 0" 4 (A r ), one traces forward for Q A (Np), and 
then one can easily solve the set of ODEs for Q A (N) via using the initial 
conditions Sl A (N F ) = N^ B & B (N F ). 



Substituting the various numerical solutions for all the relevant quantities in (]3! 
one obtains the values for all the three non-linearity parameter f^i, tnl and gNL- 



5.2 Numerical Results 

Now, we present the numerical results applying the strategy discussed so far 
for our potential f l25]) . The various non-linearity parameters /nl,tnl and gNL 
have been estimated for slow-roll as well as beyond slow-roll regime for the four 
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trajectories mentioned before corresponding to sufficiently large e-folding N > 50. 
We observe that in the slow-roll regime, the /jvl (as well as r NL ,g NL ) are very 
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Table 3: Various non-linearity parameters estimated up to the end of slow-roll 
region. 

small which is also something quite expected [52]. Table [3] presents the collection 
of non-linearity parameter values estimated in the slow-roll regime. 
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Table 4: Various non-linearity parameters estimated in the beyond slow-roll 
regime. The time corresponding to A 7 "' is simply the time where axion p w first 
crosses its minimum and oscillates up to the end of inflation. 

Table H] shows that one can get large non-linearity parameters in the beyond 
slow-roll region for these trajectories except the single- field one. Especially, most 
of these large values are generated after the time when axion crosses its mini- 
mum for the first time, i,e. between Np and . One can see that in the beyond 
slow-roll regime (Np — N = 2 or at most 3), the trajectories acquire quick turns 
via oscillations in axionic directions when one approaches towards the end of 
inflation. This might be a reason for these non linearity parameters getting en- 
hanced to significantly large values. Another reason for having large values of 
non-linearity parameters could be attributed to one of the slow-roll parameters 
getting significantly large, rj ~ (9(10 2 — 10 3 ). Moreover, contributions to the e 
parameter com ing from each of the two-fields are such that €\ ^> e 2 in beyond 
slow-roll regima 10 l. These can cause large enhancements in the intermediate quan- 
tities such as G" 4 in the regime where the slow-roll condition is strongly violated. 
The reason could be simply thought of to be as Q A (N) = A A B (N, N')N B (N') and 

10 In two-field model with t\ 3> £2, it has been argued that Jnl can be enhanced by a factor 
of 0(ei/e 2 ) to its naively expected 0(e) value Note that such hierarchies in slow-roll 

parameters can be extremely crucial and this has been a key in realizing large /jvl even in 
slow- roll regime [S3]- Similar example could be a two-field DBI inflationary model like [55]. 
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Ag(JV, iV') are defined through exponential of integrals involving Pjg[N) which 
depends on the combinations of slow-roll parameters. However, in the slow-roll 
regime, A^ B (N, N') ~ 0{1) and this results in the /jvl parameters to be of slow- 
roll suppressed values for various trajectories. Also, in [5H], it has been argued 
that the value of Jnl parameter can be expected to be crucially enhanced or 
suppressed if A A B (N, N') happens to be so. 

6 Conclusions and Discussions 

In this article, we generalized the standard single field poly-instanton inflationary 
model to a two-field model in which inflation is driven by a combined dynamics 
of a (Wilson) divisor volume mode and the respective C 4 axion. In this setup, we 
have focused on two aspects. The first one was the study related to the two-field 
inflationary model. The second one has been an investigation of the possibility 
to realize the non-Gaussianities by computing the non-linearity parameters such 
as JmliTnl and gNL- This was done in the slow-roll as well as in the beyond 
slow-roll regime. 

In the context of the first aspect, we studied the evolutions of background 
fields and explored the various possible types of inflationary trajectories. We 
observed that, depending on the choice of initial conditions, one can have infla- 
tionary trajectories corresponding to a wide range (order one e-folding to quite 
large) number of e-foldings. However, we mainly focused on the class of trajec- 
tories which could produce order 50 (or more) e-foldings. The nature of various 
trajectories are also quite different; some are attracted (repelled) inwards (out- 
wards) to the attractor point in a straight line and hence implying that such 
trajectories have no isocurvature perturbations. The other type of trajectories 
have significant curving due to the axion dynamics and do have isocurvature 
perturbations. Further, if the axion initial conditions are not very far from the 
respective minimum, the respective trajectories are attracted into the nearest 
valley. Interestingly, there are some trajectories which are predominantly axionic 
before getting trapped into an attractor point. This situation is similar to the 
rotation of a roulette ball before being trapped into a particular slot, and so jus- 
tifies the name "roulette" inflation [JB]. We analyzed the running of the slow-roll 
parameters, the Hubble rate and the power spectrum of scalar perturbations in 
terms of e-foldings. Similar to the single-field inflationary models [281 E2] realized 
in type IIB LARGE volume orientifold compactifications, the values of slow-roll 
parameters e and rj are quite hierarchial. The slow-roll parameter e is very small 
(e ~ 10~ 9 ) during slow-roll and it increases only up to order 10 -6 values when 
slow-roll conditions are violated by r\ ~ 0(1). In fact, the e parameter increases 
exponentially fast near the end of inflation and at the same time the r\ param- 
eter also gets enhanced to a very large values. Thus, from the point of view of 
e-foldings, a very narrow window is available in the 'beyond' slow-roll regime in 
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which rj parameter is very large. Moreover, a significant hierarchy ei ^> €2 is also 
realized within the components of e in beyond slow-roll regime. Also, as we have 
found several trajectories with significant curving, it has been a good motivation 
for looking at the non-Gaussianities signatures in this "roulette poly-instanton 
inflation" setup. 

Subsequently, we have investigated the possibilities for generating detectable 
values for the non-linearity parameters /jvl, r NL and g^L for various trajectories. 
We followed the strategy developed in [HSl EH] which is applicable for a given 
generic scalar potential and is also valid in the beyond slow-roll regime. For 
the trajectory which has initial conditions such that axion sits at its minimum 
position to start with, this two-field inflationary process reduces to its single- 
field analogue [32]. Thus, for this trajectory, we get quite small value for the 
non-linearity parameters /jvl, t nl and <?jvl in the slow-roll regime. In fact, 
within the slow-roll regime, we obtain small values for each of the the three non- 
linearity parameters for all the 'types' of trajectories we have classified in this 
setup. However, in the beyond slow-roll regime (Np — N = 2 or at most 3), 
the trajectories acquire turns due to oscillations in axionic directions when one 
approaches towards the end of inflation. These might be responsible for the non 
linearity parameters getting enhanced to significantly large values. Moreover, a 
sharp increment of the rj parameter towards the end of inflation could be another 
reason for the non-linearity parameters getting enhanced in the beyond slow-roll 
regime. Recently, in the context of estimating non-Gaussianities signatures in the 
beyond slow-roll regime, large /jvl values have been observed in [22 ED] with a 
potential having a combination of various exponential terms such that the Hubble 
rate is sum separable. Such a potential has been argued to be possibly realized 
in string models such as [28]. Although, our two-field scalar potential does not 
have a separable form for the Hubble rate (and hence one can not expect to 
reproduce our results from the strategy developed in [591(60]), nevertheless there 
are important similarities which could be responsible for the large values of /jvl 
parameter. A few of these are, 

• Large /jvl is realized only in the region of the final stage of inflation where 
there is no significant increment in the number of e-folds. 

• In both of the cases, there is a significant enhancement in the //-parameter 
towards the end of the inflation resulting in hierarchial values for the e and 
77 parameters. 

• Without going into all the details of [60], we state as an observation that 
the exponent appearing in the form of the Hubble rate as e~ ak< ^ k results 
in /jvl ~ 0{nak/ y/rnk)-, where n is the number of fields involved in the 
dynamics and o^, are some model dependent parameters. For ~ 
0(100), nik ~ 1, one gets /jvl ~ 0(100). Now for our case, after considering 
the canonically normalized forms of the divisor volume modulus, as seen in 
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single- field analogue in [32] , the exponential term in the potential appears 

4/3 2/3 

as e _Qfe fc where ~ V and V is the stabilized value for the Calabi 

Yau volume. For the present case, V ~ 905. Moreover, recall that such 

_ ,4/3 

exponential form results in r] ~ ^ ~ e~ afc ^ which is similar to those of 

M- 

Along these lines, our analysis could be a supporting step towards looking for 
the non-Gaussianities in the beyond slow-roll limit. Also, a systematic study of 
what happens after the end of inflation, for example following the motion/ decay of 
inflaton fields and reheating issues etc. can be extremely crucial in such modelJ^l. 
In the context of Kahler moduli inflation models [281 EQ], it has been found 
that there is a possibility of inflaton dumping all its energy into the hidden 
sector instead of the visible sector [HH E2]- For investigating this aspect, one 
has to extend our setup by explicitly embedding the visible sector. It would be 
interesting to know how generic is the problem and we hope to get back with this 
aspect in a future work. 
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A Collection of the relevant expressions 

Recalling the field redefinitions, 

the background field-evolution is governed by the following equivalent expressions, 

m <Pi = Fi ^2 = F 2 = -3^ 2 - — 

where D is the covariant derivative defined as Dtp?; = d(p2+Y a i )C ip b 2 dip1 1 and Hubble 
rate is defined via H 2 = | (V + | Qabip^ip^) ■ Now, utilizing the aforementioned 



n We thank E. Copeland and A. Mazumdar for pointing this out to us and bringing our 
attention to the related studies in [501 13T1 1521 155] 
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relations, the various useful expressions appearing at the intermediate stages in 
the computation of non-linearity parameters (/jvio t~nl and Qnl) can be easily 
derived. The expressions for Pg = (j^jrj , „ are simply given as 



at y*=vf 0) {N) 



P a 2 h = + 777^3 VaV * - ^cM f2 ^2 (41) 

pa ib = jj$b ~ ^jp (Gbdvt) 

The expressions for Q^ c = (~^Ar J are given as 

Q a \l = -^3 ^ H c + ^ ^ V b V c - if w ^ 

Q a ll = ~ 6a c V b + 12^B ^ Wcdf^) 

-^- 3 S a c (G bd ^)-^^ a 2 G bc (42) 
V a V b V c - V c R a pbq tftp* 



12 H 5 



Q a Z = ^ V a b (Gdtfi) - ^ V a V b (Gafl) - 2R\ d <p d 2 

Q a lL = ^ K (gupi) - y^JjE ya (Ml) K - R a dcb <£ 
Q a ltc = -j^jp v a (G b M) (G b M) + ^v a g bc 

After including the slow-roll corrections in A AB [ST] and keeping in mind that 
background metric for our case is diagonal, we get the following components, 



l + 2e + «l -4e + £U ^- ) - g bb [ ^- ) -2 G 



Aaa pa 

n -y |- , ~ , - ^ - , ~™y dN j ■>«. y d 

A°\ = for a ^ b, where a = {1, 2} and a ~ 0.7 
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,«6 Q ac v c Q hd v d G ac G bd (Vcd - v q cd v q ) , ab 

n i2 — TP) 77 ~~ A 2i) l 4t) J 



^22 



V 2 V 

7° V c G ak (V ck - T\ k V q ) \( V c V b G cp G bq (v pq - T^Vi) 



V 2 V \ V 2 V 



In the first two-lines of (|43[) . the indices are not summed over. Furthermore, the 
expressions for H^, Hj£, H^bc ar e elaborated as, 

Hi = ^ V a , Hi = (g ab cp b 2 ) (44) 
H * = 6H Vab ~36lp VaVb 

H ab = ~^Jp Va (Gbc<P2) 
H ab = - gg jj 3 Vb {Gac^D 

H °- b = 6# ^ ab ~ 36 H 3 (GacV'di.Gbdyi) 

H^ c = g^r V abc - {V ab V c + VacVt + KHc) + y^Jp V a V b V c 

H lbc = ~^J[3 V ab(Gcd<P 2 ) + f^Jp VaVbiGcdVi) 
H lbl = VaiGbdf^Vc ~ 7^~Fp VaciGbd^t) 

H abc = V a{Gbd<pt)(Gcp¥%) - ^Jp Va Gbc 



H lbl = Y^Jp (Gad^) VbV c ~ ^Jjp (Gad^) Vbc 



H abc = ^Jp (Gad^iWb (GcpVl) ~ ^Jp V b 

H abc = ~^Jp Gab V c + J^Tp (GadVt) (Gbp^l) V c 

H abc = J2Jp (GadVt)(Gbp¥ P 2 ) (GcqVl) ~ ^Jp Gab (Scdft) 

■ Gac (Gbd^i) — 7Tn Tr i (Gad^i) Gbc 



The expressions of various derivatives of e-folding N evaluated at a final time- 
hypersurface tp (e.g. N%, N^ B , N% BC ) which is used for providing the initial 
conditions while solving for the ODEs backward in time are given as, 



H V F?> 



at y=ip(°){N F ) 
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AB 



U M 



N 



ABC 



Zabc 
H V FV 



at (p=<p(°)(N F ) 



where 



Uab = H ab + 2 (H c P C A + F c H CA ) Ni 



+ (F C H CV F v + H c P% F v ) N F A N F 



Zabc = H AB c + 



H v (Q V ZT F £ + Pf P £ ) F T + Hvejr F V F £ F T 



+3F V H V£ P*F 



N A NgNc + 3 



{H AV£ F v + H AV Pf) F £ 



+2 F V H V£ P £ A + H v (Qf A F £ + Pf P £ A ) 
+3 {2H AV P^ + F v H VAB + H v Q V AB ) N* 



N§Nf 



+3 {F v H V£ F £ + H v Pf F £ ) N F N F C + 3 {F v H va + H v P%) N x 
Finally, the expressions for Q% cv = 

1 



BC 



at ^=<pf Q) (N) 



Q 



alll 
lbcd 



12H 5 



V2 He (Gdptfl 



1 



6H 3 



s a d v t 



be 



72H 7 



Q 



oll2 



6# 3 



d Vbc 



72H 7 



<f a 2 v b v c (g dp <4) 

5 a d V b V c 



12H 5 



v a 2 v b v c (g dp <p p 2 ) 



12H 5 



12H 5 



S a d V h V c 



<^2 V bc (Gdp<f%) 



Q 



al21 
lbcd 



S\V hd + -^—5\V b V d 



6H 3 



12H 5 

5 1 



Q 



o211 



72# 7 

5 °& ^ c Vd ~ 6^ ^ Vcd 

1 5 

+ ^2 (Gbp^fl) Vd ~ ¥>2 (GbpVfi) V c V d 



12H 5 



f^a\22 



5 1 



29 



12H 5 r2 ^ oc " a ' 12H 5 u b p ^ 2 ^ Vd 



&\ v c {g dP rt) + -Lr 5 a d {g^d v c 



12H 5 



12H 5 
1 



12H 5 

Pi (Gbp<f%) V c {Gdqtpl) 



1 



P 2 Vb Gcd 



6H 3 



72H 7 ^ zy p ^ z ' i2i/5 

<*°6 (Sc^f) (Sd^l) + S\ (Gdqtpl) + V?2 Gic {Gdqtpl) 

5 



5 a b Gcd + 5 a c Gbd + 5 a d Gbc 
1 



72H 7 



P 2 (Gbp^pDiGcr^DiGds^Pl) 



12H 5 



^diGbqtpDiGcrtPz) + ^2 G b d(Gcr-tp 2 ) + ^2 (.Gbrfl) Gcd 
1 



feed 



6# 3 



^ a fec v d + V° M V c + V\ Vcd + v™ cd H + v\ v u 



-v a d v bc + v a v bcd 



1 



12# 5 

- v° He v d + v a v bd v c + ^ a v b v cd 



v\ v c v d + v a c v b v d + v a d v b v c 

5 



72H 7 



v a v b v c v d 



1 n c (M)-^n^M) 



6# 3 



QH 3 



12H 5 



V a c V b {GdM) 



12H 5 

V\ d (Gc^l) 



V a V bc (G dq <f q 2 ) + V a V b V (G dq <p q 2 ) 



72H 7 

1 V\ {Gcqtfl) V d+^jp Va V b (Gcqfl) V d 



12H 5 



12H 5 



V 1 (Gbqtfl) V cd 



+ 



V a (G bq <p 9 2 )V c V d 



n 19wp\) Vc 



72H 7 ^"iru c a 12R5 
^3 V\ Gcd - ^ V\ (Gcp<{%) (Gdq^l) 

V a V b Gcd + V a V b (G CP P P 2 ) (Gdqtpl) 



12H 5 



72H 7 



V a Gm V - 777777 V a c (G bp p p 2 ) (GdqPl) 



12H 5 



72H 7 



12H 5 

V a (G bp p p 2 )V c (GdqP q 2 ) 



QH 3 



V a r G 



Ixl 



- n (G bp tp p 2 ) (Gdqtpl) + ^V a (G bp tp p 2 ) V (Gdqtpl) 



12H 5 



QH 3 



v a c G 



bd 



12H 5 
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72H 7 
V a G bd V c 



K 6U 1 

^ V a {G bp p P 2 ) (GcqPl) (GdrPl) ~ J^r V° ' G b d (Gc P ^ 2 ) 



(47) 



12H 5 



V a (G bp tp p 2 ) Gcd 



1 



12# 5 



V a G bc (Gd P tp\) 
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